Comparison of transcriptome responses of the liver, tail fin, and olfactory epithelium of Rana [Lithobates] catesbeiana tadpoles disrupted by thyroid hormones and estrogen.

Thyroid hormones (THs) are important developmental regulators in vertebrates, including during the metamorphosis of a tadpole into a frog. Metamorphosis is a post-embryonic developmental period initiated by TH production in the tadpole thyroid gland. The two main bioactive forms of TH are L-thyroxine (T4) and 3,5,3'-triiodothyronine (T3); these hormones have overlapping but distinct mechanisms of action. Premetamorphic tadpoles are highly responsive to TH and can be induced to metamorphose through exogenous TH exposure, making them an important model for both the study of vertebrate TH signaling and endocrine disrupting chemicals (EDCs). It is important to differentiate TH-mediated responses from estrogenic responses in premetamorphic tadpoles when assessing dysregulation by EDCs as crosstalk between the two endocrine systems is well-documented. Herein, we compare the RNA-sequencing-derived transcriptomic profiles of three TH-responsive tissues (liver, olfactory epithelium, and tail fin) in premetamorphic bullfrog (Rana [Lithobates] catesbeiana) tadpoles exposed to T3, T4, and estradiol (E2). These profiles were generated using the latest available genome assembly for the species. The data indicate that there is a clear distinction, and little overlap, between the transcriptomic responses elicited by E2 and the THs. In contrast, within the THs, the T3- and T4-induced transcriptomic profiles generally show considerable overlap; however, the degree of overlap is highly tissue-dependent, illustrating the importance of distinguishing the two THs and the affected signaling pathways within the target tissue type when evaluating hormone active agents. The data herein also show that E2 and TH treatment can uniquely induce significant changes in expression of their respective "classic" bioindicator transcripts vtg (E2) and thra, thrb, and thibz (THs). However, care must be taken in the interpretation of increased vep or esr1 transcripts as a change in transcript levels can be induced by THs rather than solely E2.

Relationship between serum thyroid hormones and their associated metabolites, and gene expression bioindicators in the back skin of Rana [Lithobates] catesbeiana tadpoles and frogs during metamorphosis.

Anuran metamorphosis is characterized by profound morphological changes including remodeling of tissues and organs. This transition is initiated by thyroid hormones (THs). However, the current knowledge of changing levels of THs during metamorphosis relies on pooled samples using methods known for high variability with sparse reporting of measured variation. Moreover, establishing a clear linkage between key gene expression bioindicators and TH levels throughout the metamorphic process is needed. Using state-of-the-art ultra-high performance liquid chromatography isotope-dilution tandem mass spectrometry, we targeted 12 THs and metabolites in the serum of Rana [Lithobates] catesbeiana (n=5-10) across seven distinct postembryonic stages beginning with premetamorphic tadpoles (Gosner stage 31-33) and continuing through metamorphosis to a juvenile frog (Gosner stage 46). TH levels were related to TH-relevant gene transcripts (thra, thrb, and thibz) in back skin of the same individual animals. Significant increases from basal levels were observed for thyroxine (T4) and 3,3',5-triiodothyronine (T3) at Gosner stage 41, reaching maximal levels at Gosner stage 44 (28 ± 10 and 2.3 ± 0.5 ng/mL, respectively), and decreasing to basal levels in juvenile frogs. In contrast, 3,5-diiodothyronine (T2) increased significantly at Gosner stage 40 and was maintained elevated until stage 44. While thra transcript levels remained constant and then decreased at the end of metamorphic climax, thrb and thibz were induced to maximal levels at Gosner stage 41, followed by a decrease to basal levels in the froglet. This exemplifies the exquisite timing of events during metamorphosis as classic early response genes are transcribed in anticipation of peak TH concentrations. The distinct T2 concentration profile suggests a biological role of this biomolecule in anuran postembryonic development and an additional aspect that may be a target of anthropogenic chemicals that can disrupt anuran metamorphosis and TH signalling. Hence, as a second aim of the study, we set out to find additional bioindicators of metamorphosis, which can aid future investigations of developmental disruption. Using a sensitive nanoLC-Orbitrap system an untargeted analysis workflow was applied. Among 6,062 endogenous metabolites, 421 showed metamorphosis-dependent concentration dynamics. These potential bioindicators included several carnitines, prostaglandins and some steroid hormones.

Antimicrobial peptides from Rana [Lithobates] catesbeiana: Gene structure and bioinformatic identification of novel forms from tadpoles.

Antimicrobial peptides (AMPs) exhibit broad-spectrum antimicrobial activity, and have promise as new therapeutic agents. While the adult North American bullfrog (Rana [Lithobates] catesbeiana) is a prolific source of high-potency AMPs, the aquatic tadpole represents a relatively untapped source for new AMP discovery. The recent publication of the bullfrog genome and transcriptomic resources provides an opportune bridge between known AMPs and bioinformatics-based AMP discovery. The objective of the present study was to identify novel AMPs with therapeutic potential using a combined bioinformatics and wet lab-based approach. In the present study, we identified seven novel AMP precursor-encoding transcripts expressed in the tadpole. Comparison of their amino acid sequences with known AMPs revealed evidence of mature peptide sequence conservation with variation in the prepro sequence. Two mature peptide sequences were unique and demonstrated bacteriostatic and bactericidal activity against Mycobacteria but not Gram-negative or Gram-positive bacteria. Nine known and seven novel AMP-encoding transcripts were detected in premetamorphic tadpole back skin, olfactory epithelium, liver, and/or tail fin. Treatment of tadpoles with 10 nM 3,5,3'-triiodothyronine for 48 h did not affect transcript abundance in the back skin, and had limited impact on these transcripts in the other three tissues. Gene mapping revealed considerable diversity in size (1.6-15 kbp) and exon number (one to four) of AMP-encoding genes with clear evidence of alternative splicing leading to both prepro and mature amino acid sequence diversity. These findings verify the accuracy and utility of the bullfrog genome assembly, and set a firm foundation for bioinformatics-based AMP discovery.